Character of magnetic excitations in a quasi-one-dimensional antiferromagnet near 
the quantum critical points: Impact on magneto-acoustic properties 
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We report results of magneto-acoustic studies in the quantum spin-chain magnet MCI2- 
4SC(NH2)2 (DTN) having a field-induced ordered antiferromagnetic (AF) phase. In the vicinity 
of the quantum critical points (QCPs) the acoustic C33 mode manifests a pronounced softening 
accompanied by energy dissipation of the sound wave. The acoustic anomalies are traced up to 
T > Tn, where the thermodynamic properties are determined by fermionic magnetic excitations, 
the "hallmark" of one-dimensional (ID) spin chains. On the other hand, as established in earlier 
studies, the AF phase in DTN is governed by bosonic magnetic excitations. Our results suggest the 
presence of a crossover from a ID fermionic to a 3D bosonic character of the magnetic excitations 
in DTN in the vicinity of the QCPs. 

PACS numbers: 72.55.+S, 75.45.+J 



The interest in quasi- ID quantum spin systems has 
grown considerably during the last decade. This is fos- 
tered by the progress in preparing materials with well- 
defined ID spin subsystems and the possibility of an- 
alyzing the experimental data with the help of non- 
perturbative theories for ID models.- In addition, such 
systems often manifest quantum phase transitions at 
T=0 which are governed by parameters other than the 
temperature. True ID models do not exhibit any long- 
range order at finite temperatures.- Real quasi-lD anti- 
ferromagnetic (AF) materials, containing weakly coupled 
spin chains with gapless spectra of their low-lying exci- 
tations, are usually magnetically ordered at low temper- 
atures. At temperatures higher than the Neel tempera- 
ture, Tn, but of the order of the exchange constant, these 
systems behave as quantum spin chains, where any long- 
range magnetic order is destroyed by enhanced quantum 
fluctuations^ One should note that quasi-lD magnets, in 
which the low-energy eigenstates of their ID subsystems 
have spin gaps, usually do not manifest long-range mag- 
netic ordering. 2 However, an external magnetic field can 
close the spin gap, A, and for H > H c ~ Aa quantum 
phase transition to a phase with gapless spin excitations 
takes place. A further increase of the field yields a sec- 
ond quantum phase transition to a spin-polarized phase 
at H > H s . In the spin-polarized phase the low-energy 
excitations are also gapped. Hence, the magnetically or- 
dered phase can be observed in the field domain where 
spin excitations are gapless, and the Neel temperature in 
such systems is field dependent. The magnetic suscepti- 
bilities of a quasi-lD spin system in mean-field approxi- 
mation can be written as 



(Xq) (1) (T) 



l-ZJ ± (q)(v«)(D(T) 



(1) 



where the superscript (1) denotes the susceptibility of 
one chain, a = x,y, z, J± is the weak interchain exchange 
constant, Z is the coordination number, and q is the wave 
vector. The quasi- ID spin system becomes ordered when 
the denominator becomes zero (which defines Tn{H)). 

The low-T thermodynamics of a state with long-range 
magnetic order is determined by bosonic excitations, 
magnons. Recently, several groups have observed phe- 
nomena in some AF systems that have been interpreted 
as Bose-Einstein condensation (BEC) of magnons, viz., 
as a thermodynamically large number of magnons in the 
same ground state^i^ For quasi-lD spin systems at 
T > Tn the projection of a single spin may have only 
a limited number of values (e.g., two values for spin-1/2 
systems, three values for spin-1, etc.). That is why ther- 
modynamic properties of, e.g., AF spin-1/2 chains are 
often determined by low-energy eigenstates which be- 
have as interacting fermionsX The fermionic nature of 
the low-energy excitations of these quantum spin chains 
with only short-range correlations is related to the lim- 
ited number of projections of each spin. Very recently 
it has been shown£ that the behavior of a spin-1 spin- 
gapped system in magnetic fields close to field-induced 
quantum critical points (QCP's) can be described by free 
fcrmions as well. The fermionic nature of these excita- 
tions is also related to the limited number of projection 
values at each spin-1 site. The fermionic behavior of ex- 
citations is characteristic for quantum spin chains with 
short-range correlations. From this perspective it is very 
interesting to study the behavior of a spin-gapped quasi- 
lD AF system in the vicinity of Tn, close to H c and 
H s . Here, low-lying excitations of the quasi-lD system 
should change their statistical properties from fermionic, 
at T > T N (H), to bosonic, at T < T N (H). Hence, by 
varying H and T one may observe features in the same 
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spin system characteristic either to fermions or bosons. 

One of the best candidates for studying such a 
crossover in the excitation statistics is the spin-1 system 
dichloro-tetrakis thiourea- nickel (II), NiCl 2 -4SC(NH 2 )2, 
known as DTN. Recently, some features in the magneti- 
cally ordered phase of DTN at H > H c and T < T N (H) 
were interpreted as BEC of spin degrees of freedom.— The 
bosonic character of the spin excitations in DTN in the 
magnetically ordered phase has been corroborated and is 
considered to be a well established fact. In this work, we 
study magnetic and magneto-acoustic characteristics of 
DTN near the critical values of H c and H s . We show that 
the behavior of the observed properties outside of the 
AF phase can be well described by an effective fermionic 
model of low-lying spin excitations. In this way DTN 
manifests ID fermionic character of spin excitations at 
T > Tff(H). This fact, together with previous results, 
showing bosonic 3D behavior of magnetic excitations in 
DTN for T < T N (H)£ leads us to the conclusion that a 
crossover from fermionic to bosonic features of the low- 
lying magnetic excitations takes place at Tjv near the 
quantum critical points. 

Following Ref. @ we describe the spin-1 chain at low 
excitation densities using an effective free-fermion theory 
with two branches of low-energy states. Two branches 
are used because the strong single-ion "easy-plane" mag- 
netic anisotropy D observed in DTN 8 splits the spin 
triplet of the spin-gap modes, and makes one of them inef- 
fective at the critical fields. 6 The two fermionic branches 
have features at H c and H s , while for D <C T the contri- 
bution of the third branch is exponentially small and can 
be neglected in our approximation. Both critical fields 
are related only to the lowest branch of our model (they 
correspond to van Hove singularities, connected with two 
edges of that branch). However, the field dependence is 
present also in the temperature-dependent factors of both 
branches. In Fig. 1 the solid line shows the calculated 
field dependence of the magnetization of the quasi- ID 
spin system (ID subsystems are considered within this ef- 
fective free-fermion model) for T slightly above the phase 
boundary Tn(H) (see inset in Fig. 1), where the suscep- 
tibility of the quasi-lD system diverges. For comparison 
we also plot experimental data taken at T = 0.6 K. 9 
Note that T^(H) < 0.6 K in the vicinity of the quantum 
critical points H c and H s , since there is a line of phase 
transitions T N {H) with T N (H C ) = T N (H S ) = and the 
system in our model is not in the magnetically ordered 
phase inside of the interval H' c < H < H' s . H' c and H' s 
are the critical fields at non-zero T. There is a good 
agreement between the effective free-fermion theory and 
the experimental data at H < H' c , H > H' s , and near 
the critical values of H. On the other hand, inside of 
the interval H' c < H < H' a , the real system is ordered, 
Tn(H) > 0.6 K, and our ID fermionic description cannot 
be applied. 

Ultrasonic investigations are a powerful experimental 
technique to study various phase transitions and critical 
phenomena. This technique is well established as an im- 
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FIG. 1: (Color online) Low-temperature magnetization (T = 
0.6 K) of DTN as a function of external magnetic field 
H || [001] (circles) .— The line is the result of the free-fermion ef- 
fective theory. The inset sketches the temperature - magnetic 
field phase diagram of DTN.— The quantum critical points are 
H c ~ 2.1 T and H a ss 12.6 T. The maximum temperature of 
the AF order is T% ax « 1.2 K. 



portant tool for the investigation of low-dimensional spin 
systems.- Spin-lattice interactions are responsible for 
the attenuation of acoustic waves and influence the sound 
velocity in magnetic crystals. These interactions are con- 
nected either with a strain modulation of the exchange in- 
teractions or with a magnetostrictive coupling of a single- 
ion typeJ^ We have performed measurements of the rel- 
ative change of the sound velocity and attenuation in 
DTN, using a phase-sensitive detection technique based 
on a standard pulse-echo method with a set-up similar 
to the one described in Ref. 10]. DTN has a tetragonal 
crystallographic symmetry (space group 14) with two for- 
mula units in the unit cell. The investigated single crystal 
has a size of about 2 x 2 x 4.1 mm 3 . Since the as-grown 
surfaces of the crystal were smooth and parallel, we glued 
piezoelectric film transducers directly to the surfaces nor- 
mal to the crystallographic [001] direction, without any 
additional sample polishing. This geometry corresponds 
to the longitudinal acoustic C33 mode, with propagation 
direction and polarization along the spin chains. A num- 
ber of ultrasonic echoes have been detected. The absolute 
value of the sound velocity at liquid-helium temperature 
has been determined as vi — 2640±20 m/s. Note that the 
measurement accuracy for a relative change of sound ve- 
locity is of the order of 10~ 6 . The sample- length change 
is relatively small for the applied temperatures and mag- 
netic fields^ Therefore, we did not have to take into 
account any length-change corrections to the sound ve- 
locity. The data have been collected using the ultrasonic 
signal at 78 MHz. The magnetic field was applied along 
the [001] direction, i.e., parallel to the sound-propagation 
direction. Figure 2 shows the magnetic-field dependence 
of the relative change of the sound velocity and atten- 
uation of the C33 mode in DTN for T below the maxi- 
mum of TJ^ ax w 1.2 K. There is a pronounced softening 
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FIG. 2: (Color online) Field dependence of the relative change 
of the sound velocity (top) and of the sound attenuation (bot- 
tom) of the acoustic c 33 mode in DTN at T below T$ ax . The 
magnetic field was applied along the [001] axis. The ultrasonic 
frequency was 78 MHz. The insets show the sound velocity 
and attenuation in the vicinity of H c in enlarged scale. 



of the C33 mode in the vicinity of both critical fields, 
though the anomaly at H' s is approximately one order 
of magnitude larger than that at H' c . There is a rela- 
tive increase of Av/v = 7xl0 -4 between the sound ve- 
locity at H = and H > 12.6 T, where all spins are 
polarized. The relative decrease of the sound velocity 
reaches about 4xl(T 3 at 12 T and 0.3 K. The soften- 
ing of the C33 mode is accompanied by a peak in the 
sound attenuation. Both the sound-velocity and sound- 
attenuation anomalies become smaller and broader with 
increasing T. The H dependence of the sound velocity in 
the ordered phase (far from the critical regions) resem- 
bles c-axis magnetostriction dataiii However, the change 
in the sound velocity cannot be explained by the lattice- 
parameter change, since the length change observed in 
Ref. [Ill is too small. In Fig. 3, we show the field depen- 
dence of the sound velocity and attenuation of the C33 
mode in DTN at various temperatures above T^ ax . One 
can see some transformation of the acoustic anomalies by 
moving from T < T™ ax to T > T™ ax . Here, the soft- 
ening of the C33 mode disappears in the vicinity of H c ; 
only a smooth increase in the sound velocity is detected. 
Close to H s one can still observe a minimum in the sound 
velocity and a maximum in the sound attenuation, but 
those anomalies are smaller in amplitude and broader 
than the corresponding ones measured below T^ ax (H) 



(Fig. 2). 

In magnetic materials the dominant contribution to the 
spin-lattice interactions mostly arises from the exchange- 
striction coupling. In our calculations we assumed that in 
DTN the spatial dependence of the magnetic anisotropy 
constant is weaker than the spatial dependence of the 
exchange integrals. In this case, one can expect that 
only longitudinal sound waves interact with the spin sub- 
system. According to Ref. [lU, the relative renormaliza- 
tion of the longitudinal sound velocity can be written as 
(Av/v) = -(Ai + A 2 )/(iVwk) 2 , where 

A 1 =2|Gg(k)| 2 (^) 2 X ^+T^ J2 l^(k)| 2 (^) 2 , 

q a—x,y,z 

A 2 =H*(k){Stf + ^J2 E ff q( k )*q- ( 2 ) 

q a— x,y,z 

Here, N is the number of spins in the system, = vk is 
the low-fc dispersion relation with sound velocity v in the 
absence of spin-phonon interactions, (S§) is the average 
magnetization along the direction of the magnetic field, 
Xq V ' z are non-uniform magnetic susceptibilities, and the 
subscript corresponds to q = 0. In the framework 
of our effective free-fermion model the temperature and 
magnetic-field dependence of the uniform susceptibility 
of one spin chain can be written as 

/ Z \(i) _ f Hs 

1 + cosh(g/r) cosh(.T/r) 

[cosh(i//T) + cosh(x/T)] 2 ' (/>1 

where we set the units for the effective g-factor, Bohr's 
magneton, and Boltzmann's constant equal to 1. For 
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FIG. 4: (Color online) Attenuation (upper surface, blue) and 
relative change of the velocity (lower surface, red) of the lon- 
gitudinal sound versus H and T', calculated in the framework 
of the proposed theory (arbitrary units were used for all pa- 
rameters, see text for details). 



spin systems with AF interactions the main contribution 
to the summation over q in Eqs. (2) comes from terms 
with q — 7r, 



06) 



(1) 



Ha 



dx 



KJH C y/(H2-X*)(x* 

sinh(x/T) 



cosh(H/T) + cosh(x/T) 



(4) 



To calculate magnetic susceptibilities of the quasi-lD 
spin system we use Eqs. (1), (3) and (4). 

The renormalization is proportional to the spin- 
phonon coupling constants 
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(5) 



Here, m is the mass of the magnetic ion, J^ n denote 
(anisotropic, generally speaking) exchange integrals, e k 
is the polarization of the phonon with wave vector k, 
and R„ is the position vector of the n-th siteJ^ Figure 4 
(lower surface) shows the H and T" dependence of the rel- 
ative change of the longitudinal sound velocity of a quasi- 



lD spin system calculated in the framework of the effec- 
tive free-fermion model. We fixed H c and H s and used 
arbitrary units for T" in Fig. 4 (they are not equal to the 
temperatures in the experiment). It is challenging to cal- 
culate Xq' v m the framework of the used model. Clearly 
they have to be smooth functions of H and T, except at 
the line T N {H). G"(k) and i?£(k) in Eqs. (2) are also 
unknown for any a (one of the coupling constants can be 
estimated using Ref. 11). That is why, in order to obtain 
the results presented in Fig. 4, we used J±(q) = 0.18 from 
Ref. H, Xq=o an d Xq=-w multiplied by some (not known) 
values of the spin-phonon coupling constants, and £P(k) 
two times smaller than G 2 (k). The temperature of the 
divergence in the magnetic susceptibility of the quasi-lD 
system is generally determined by anisotropic couplings 
between the spin chains (these couplings are unknown). 
The divergence at Tn(H), which we used in our theory, 
does not depend on the direction of the order parame- 
ter. Such divergences are present in a quasi-lD model, 
when any component of the magnetic susceptibility (but 
with different phase-transition temperatures, T)v(-ff)) is 
considered. Even in this approximation our simplified 
theory reproduces the main features of the experimen- 
tally observed behavior. Our model reproduces the pro- 
nounced minimum at H' s , the almost field- independent 
behavior at H > H' s and H < H' c , the larger value of 
Av/v for H > H' s as compared to H < H' c) and the 
maximum (with Av/v > 0) in the interval between H' c 
and H' s . With increasing T' the features near the critical 
fields become weaker, the same way as it was observed 
in the experiment (cf. Fig. 2 and Fig. 3). At the phase 
boundary Tm{H) the susceptibility of the quasi-lD sys- 
tem diverges (see above), and our theory predicts very 
narrow and large peaks at the critical values of H (not 
shown). Therefore, for the sake of clarity, the curves in 
Fig. 4 are not plotted starting from T" = 0. T^ lax {H) in 
our units is T' = 0.02. Concerning the other values of % q 
(i.e., q 7^ 0, 7r) we affirm, as it was discussed above, that 
their inclusion does not affect the qualitative behaviour 
of the sound velocity and attenuation. Following Ref. Q~l 
we also calculated the attenuation coefficient for DTN, 



Aa(= Aa k ) 



1 

Nv 



To 



(7«?) 2 



- T E E l G q( k )| 2 (Xq) 2 



27^ 



(6) 



where 7q are the relaxation rates, which can be approxi- 
mated by 7^ = -B/TVq, where B is a material-dependent 
constant (see Ref. Il2r ). In our calculations we used the 
approximation, in which the relaxation rates do not de- 
pend on the direction and on the wave vector. The results 
are also presented in Fig. 4 (upper surface). Here, our 
theory reproduces also the main features observed in the 
experimental data: an abrupt increase of the sound at- 
tenuation near the saturation field H s , and damping with 
increasing T' . All these findings demonstrate the impor- 
tant role the fermionic magnetic excitations play in the 
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vicinity of the QCPs in DTN. 

We also tried to reproduce the observed experimen- 
tal results using the scaling-like procedure, proposed in 
Ref. Il3l . In the framework of that approach we can use 
(x%' y ) 1 j which seems more accurate than the use of 
(xl)*- 1 -* only. However, the agreement between the the- 
ory and experiment was worse than for our effective free- 
fermion model. 

Generally speaking, one could as well use a bosonic, 
say the Holstcin-Primakoff, representation of spin 
operators^ for any temperatures in DTN. However, to 
describe the behavior of spins for T > Tn(H) one has to 
take into account all interactions between these bosons 
(because the interactions are of the same magnitude as 
the energy of the free bosons), which is impossible so far. 
We do not know any other theory, bosonic or fermionic, 
which can describe the behavior of the magnetization, 
sound velocity, and attenuation in DTN better than the 
theory presented here. The situation, e.g., in weakly 
coupled spin ladders, is very different from the one in 
DTN, because in our case one cannot consider any of the 
spin-spin interactions as weak. Also, the use of hard- 
core bosons for the description of the behavior of DTN 
for T > Tm{H) cannot help because from the view- 
point of their collective behavior they can be regarded 
as fermions (i.e., only one fermion, or hard-core boson 
can be in one state). For T < T N (H) in DTN, we def- 
initely cannot use hard-core bosons (e.g., for hard-core 
bosons BEC is impossible, however, see Ref.Q for DTN). 



The advantage of our fermionic description of DTN for 
T > Tff(H) (and the mentioned bosonic description for 
T < Tn(H)), compared to the use of only a bosonic de- 
scription of low-energy spin excitations, is that in our 
approach both fermions and bosons are basically non- 
interacting. Hence, they have all features of standard 
fermions and bosons. For strongly interacting bosons, 
which is the case for hard-core bosons or the Holstein- 
Primakoff representation for T > Tn(H), one cannot, 
strictly speaking, use directly the bosonic character of 
these excitations. We finally want to note, that in our cal- 
culations we never used the symmetry of the wave func- 
tion, the other difference between fermions and bosons. 

In summary, our magnetic and magneto-acoustic stud- 
ies of the quantum spin-chain magnet NiCl2-4SC(NH2)2 
show that the behavior of the observed properties at T > 
Tn(H) can be well described by an effective ID fermionic 
model of low-lying spin excitations. This fact, together 
with previous results showing the bosonic 3D behavior 
of the magnetic excitations in DTN for T < Tjy(H)r- 
suggests the presence of a crossover from a fermionic to 
a bosonic character of the magnetic excitations close to 
the quantum critical points. The fermionic and bosonic 
nature of the magnetic excitations is related to the short- 
range correlations in the spin chains and to the long- 
range three-dimensional order, respectively. 
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